Composite Fermions in Negative Effective Magnetic Field: A Monte-Carlo Study 
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The method of Jain and Kamilla [PRB 55, R4895 (1997)] allows numerical generation of composite 
fermion trial wavefunctions for large numbers of electrons in high magnetic fields at filling fractions 
of the form v — p/(2mp + 1) with m and p positive integers. In the current paper we generalize 
this method to the case where the composite fermions are in an effective (mean) field with opposite 
sign from the actual physical field, i.e. when p is negative. We examine both the ground state 
energies and the low energy neutral excitation spectra of these states. Using particle-hole symmetry 
we can confirm the correctness of our method by comparing results for the series m = 1 with p > 
(previously calculated by others) to our results for the conjugate series m = 1 with p < 0. Finally, 
we present similar results for ground state energies and low energy neutral excitations for the states 
with m = 2 and p < which were not previously addressable, comparing our results to the m = 1 
case and the p > 0, m = 2 cases. 



I. INTRODUCTION 

The composite fermion approach pj has had a great 
number of extremely impressive successes in describing 
the physics of electrons in high magnetic fields. In 
this picture, fractional quantum Hall systems in total 
magnetic field B are described in terms of noninteract- 
ing "composite fermions" in an effective magnetic field 
B c ff — B — 2m4>on, where (f>o = hc/e is the flux quantum, 
n is the electron density, and m is a positive integer. This 
maps, for example, fractional quantum Hall states at fill- 
ing fractions of the form v = u^q/B = p/(2mp + 1) to 
integer quantum Hall states for the composite fermions 
at filling fraction i/ c s = ncj) /\B c ff\ = \p\. We will denote 
such composite fermions with 2m flux quanta attached 
to them as 2m CF. 

Jain's original approach to composite fermions 2] con- 
structed highly accurate trial wavefunctions by taking 
simple wavefunctions for the noninteracting (composite) 
fermions in the effective magnetic field, multiplying by 
Jastrow factors, and then projecting the result into the 
lowest Landau level. The early successes of this method 
were impressive 2] , despite the fact that the method was 
limited by the extreme numerical difficulty of perform- 
ing projections for systems with more than roughly 10 
electrons. 

A major theoretical breakthrough came when Jain and 
Kamilla 3] discovered a new way of writing composite 
fermion trial wavefunctions (described below), which in- 
volves a very minor modification of the projection. These 
new trial states seemed to be just as good as the originally 
proposed wavefunctions and could be numerically gener- 
ated even for systems with many electrons (40 electrons 
or more). Since that time, many important studies have 
been achieved using this method |4j, |a, LjJ. LlJ ■ However, so 
far this method has been restricted to cases where the 
effective magnetic field has the same sign as the exter- 
nal magnetic field. Results using this method have been 
published for filling fractions of the form v = p/(2mp+l) 
with p > but not for p < 0. In the current paper, we 



extend the work of Jain and Kamilla 3J so that we are 
able to handle states with p < 0. The p < states take 
more computational resources than the case of p > 0, and 
the difference in the computational resources between the 
two cases increases with the absolute value of the effec- 
tive flux. However, the computational problems turn out 
to be more severe for the smallest \p\, where the num- 
ber of particles in the system increases most slowly with 
each flux added. Fortunately, we probably need not go 
to particularly large systems to understand the physics 
of small \p\. For large \p\, describing the approach to the 
Fermi liquid-like composite fermion state, the system size 
is already large for relatively small effective flux, i.e. the 
differences of the computational requirements in the case 
of p > and p < become relatively less important. 

Although the case of negative p has not previously been 
studied for large systems, we point out that for m = 
1 the series of states with negative p and the series of 
states with positive p are essentially equivalent due to 
an exact particle-hole symmetry in the lowest Landau 
level. In fact, below, we exploit this symmetry to check 
the validity of our method. Once we have verified the 
method, we can study the properties of the m = 2 series 
for negative p and compare the results to those of the 
positive p members of this same series as well as to those 
of the m = 1 series. 

The outline of this paper is as follows. In sectionlTTlwe 
briefly review the Jain-Kamilla method. As mentioned 
above, the method has only been used previously for the 
case of positive p. In appendix lAl we show in detail how 
the crucial projection scheme of Jain-Kamilla can easily 
be generalized to handle negative p also. It is easy to see 
from the result how much additional numerical complex- 
ity is involved for negative p. In section IIIII we test our 
approach by examining the 2 CF series is = p/(2p+l). In 
particular, we examine ground state energies, excitation 
spectra, and energy gaps. We pay particular attention 
to the mass of the composite fermion and the scaling of 
the gap with p. Foro > 0, results are already available 
in the literature 0, 0, H, IE 0- F° r p < we use our 
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generalization of the Jain-Kamilla method to calculate 
energies directly, and we compare these energies to en- 
ergies obtained by particle-hole conjugating the p > 
series. Appendix 151 describes the particle-hole conjuga- 
tion transformation in depth. This comparison estab- 
lishes the accuracy of our method. In section Hvl we move 
on to examine the 4 CF series of states v — p/{^p + 1). 
Again, for p > some results are already available in 
the literature H 0, H B 0. However, for p < our re- 
sults are new. Again, we examine ground state energies, 
excitation spectra, energy gaps, and composite fermion 
effective masses. (By using particle-hole conjugation, we 
could give results for filling fractions v = 1 — p/ (4p +1)) 
similarly). We are able to make some comparison of our 
energy gaps to the experimental work of Pan et al. 8\. 

Throughout this paper we assume complete spin- 
polarization of the electrons. This should be a reason- 
able assumption for real experiments at sufficiently high 
magnetic fields. 



II. THE JAIN-KAMILLA METHOD 

Jain's original proposal^ was to construct trial wave- 
functions for fractional quantum Hall states by writing 

*jain=P{det[V' i (rS)]^ m } (1) 

where the determinant is a Slater determinant of non- 
interacting single fermion wavefunctions ipi in effective 
magnetic field B c r = B — 2mn0o, and V indicates pro- 
jection to the lowest Landau level. Here, $o is the wave- 
function of a completely filled Landau level, 



(2) 



i<3 



where Zj = (xj + iyj)/£o is the dimensionless complex 
coordinate on the plane, £q — ^Jhc/eB and the usual 
Gaussian factors exp(— j ^ \ z i\ 2 ) are understood to be 
included in the measure of the Hilbert space and will not 
be written explicitly for simplicity of notation. 

Choosing a set of single particle wavefunctions ^ to 
fill the p lowest effective Landau levels (i.e., such that 
det[ipi(rj)] represents the ground state of an integer quan- 
tum Hall state v = p), one obtains through Eq. ^ ex- 
tremely good trial wavefunctions for fractional quantum 
Hall states v — p/ (2mp+l). As discussed above, the pro- 
jection in Eq. ^is exceedingly hard to implement for sys- 
tems with more than roughly 10 electrons. For this rea- 
son, Jain and Kamilla 3] looked for an essentially equiva- 
lent formulation that would be computationally simpler. 
In their approach they begin by rewriting the wavefunc- 
tion as 



where 



* Jain =P{det [MfyJ?]} 



(3) 
(4) 



and then make the approximation that one can inter- 
change the order of projection and taking the determi- 
nant to obtain a new trial wavefunction 



JK 



det rpi(fj) 



<I> 



2m 



with 



Ur 3 ) = Jj m V {^{r 3 )J™} 



(5) 



(6) 



Although ^jk appears to be a single Slater determinant, 
it is somewhat more complicated because each ipi(fj) is 
actually a function of all of the particle positions through 
Jj . Nonetheless, this new trial wavefunction is far simpler 
to evaluate numerically. Furthermore, extensive numer- 
ical work0, 0, 0, nas shown that for small systems 
"fjK is just as good a trial state as ^jain and that both 
are extremely accurate|lfij. In the original work by Jain 
and Kamilla, it was shown how to calculate ip on a sphere 
for the case when B e ff has the same sign as the magnetic 
field B (i.e., p > 0). In Appendix lAlwe repeat the deriva- 
tion for the case where B c g has the opposite sign from 
B (i.e., p < 0). A discussion is also given there of the 
relative computational effort required to perform the rel- 
evant computations numerically. 

This technique allows one to also obtain low energy 
spectra of these fractional quantum Hall states, by simi- 
larly composite-fermionizing low energy excited states of 
noninteracting fermions as discussed in Ref. |J. 

In this paper we will perform all calculations using 
a spherical geometry^ with a monopole of charge JV^ 
flux quanta at the center. The composite fermions then 
see an effective flux 2q = Nf = N<p - 2m(N - 1). In 
the presence of this effective flux, single particle states 
are described by two quantum numbers, I and to. Here 
I = \q\ + n is the angular momentum with n — 0, 1, 2, . . . 
corresponding to the "Landau level number" or "shell" 
index, and m is the z-componcnt of angular momentum. 
A state with p filled composite fermion Landau levels 
corresponds to n = \p\ — 1. 

A low energy exciton is now formed by taking a com- 
posite fermion out of the highest occupied shell (or Lan- 
dau level) I = If = \q\ + \p\ — 1 with some mh and putting 
it in the lowest unoccupied I = + 1 shell with some m e . 
Choosing to work with states of zero total z-angular mo- 
mentum, we take the state with m e = —mh and write this 
state as \m e ). Using vector coupling (Clebsh-Gordon) 
coefficients |l2| we can construct exciton eigenstates of 
angular momentum L as|j| 



^•exciton 



= ^2 \m e )(lF,~m e :l F + l,m e \L,0) (7) 



m e — — If 



which serve as extremely accurate trial wavefunctions for 
the low energy excited states of the above discussed com- 
posite fermion ground states. 
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TABLE I: Numerical results for energies and gaps at filling fractions v = p/(2p + 1), given in units of e 2 /e£o- Calculations 
were performed using Monte-Carlo for 2 CF trial wavefunctions described in the text using 10 7 samples. Particle-hole conjugate 
pairs should give precisely the same excitation energies (and the same ground state energies once Eq. |H]is used, as it is here, 
see below). In other words, if our results were exact, the right hand columns with positive p should precisely match the left 
hand columns with negative p. Here, since we have used trial wavefunctions which are approximate (albeit extremely good), 
the agreement is not quite perfect, but it is extremely close. We note that the energies obtained by using the positive p states 
are slightly lower, which means that the trial wavefunctions with positive p (positive flux seen by the composite fermions) 
yield slightly better trial states. The number of electrons N for particle-hole conjugate pairs sums to N$ + 1, which is one 
filled Landau level. In this table, the ground state energy E g is presented for the negative p case. For the positive p case, 
in the column labeled u Eg\ P-H conj" we have put the calculated ground state energy into Eq. |H]and presented the result 
for comparison with the corresponding negative p states. The mageto-roton gap (MR-gap) is defined to be the lowest energy 
neutral excitation. The large k gap is the gap measured at the highest possible angular momentum L that we can construct 
using Eq. Q which is Z/ max = 21f + 1. All values indicated for the thermodynamic limit have been extrapolated by a simple 
linear regression over the inverse particle number using this set of data only. 



III. PARTICLE-HOLE SYMMETRY AND 
RESULTS FOR v = p/(2p + 1) 

Using particle-hole symmetry of the lowest Landau 
level, one can exactly map states at filling fraction v 
into states at filling fraction 1 — v (so long as we main- 
tain complete spin polarization). As mentioned above, 
in this paper we perform all calculations using a spheri- 
cal geometry 9]. On the sphere, the lowest Landau level 
has N,/, + 1 single particle eigenstates with the to- 
tal number of flux quanta penetrating the sphere. Thus 
states with N electrons can be precisely mapped to their 
particle-hole conjugate states with + 1 — N electrons 
(i.e. with N holes). In Appendix iBl we show that, on a 
sphere, given an eigenstate \I/ with N electrons and en- 
ergy Eqi one can write the energy Eq, of its particle-hole 
conjugate wavefunction ^ as 

Eq, = (l — / ) -Efiiicd + Ey. (8) 

where -Exiled is the energy of the completely filled Lan- 
dau level. This, of course, implies that the excitation 
spectrum of any given state is precisely the same as the 



excitation spectrum of its particle-hole conjugate state. 

We will now focus on 2 CF states of the form v = 
p/(2p + 1). The state with p is particle-hole conjugate 
of the state with p — > (— p — 1). (For example, v = 1/3 
which is p = 1 is conjugate of 2/3 which is p — —2). Ex- 
tensive numerical work has already been performed for 
positive p, calculating accurate ground state energies and 
energy gaps0, 0, IE IE 0- Using Eq. |S] this means that 
we already know the ground state energies and energy 
gaps for negative p. Here, however, we calculate these 
quantities directly using our negative p trial wavefunc- 
tions and compare to the particle hole conjugated results 
to establish the validity of our approach. 

Table [I] summarizes the numerical results for the 
ground state energies of states in this series calculated 
using Monte-Carlo|l|l|l|ll3 [13 • 

The ground-state energies we obtain using negative p 
trial wavefunctions show an outstanding agreement with 
the values obtained by particle hole conjugating positive 
p trial wavefunctions. We observed slight differences on 
the fourth significant digit, that show that the trial state 
with composite fermions in positive effective flux is very 
slightly better than the one with negative effective flux 
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FIG. 1: Quantum Hall states that are particle-hole conjugate 
of each other, are expected to have the same excitation spec- 
trum. The low energy spectra shown here for for v < i have 
been calculated with the method of Jain and Kamilla^. The 
spectra for the particle-hole conjugate states with v > i are 
calculated using the new negative effective flux trial wave- 
functions described in this paper. If the calculations were 
exact (rather than just approximate) the corresponding spec- 
tra of the particle-hole conjugate pairs would match exactly. 
Although our calculations, being based on approximate trial 
wavefunctions for 2 CFs, are not exact, we still see remarkably 
good agreement, suggesting that our new trial wavefunctions 
are approximately as accurate as the previously described trial 
wavefunctions for v < ^. The spectra shown here correspond 
(as marked) to particle hole conjugate pairs at filling fractions, 
from top to bottom 1/3 <-> 2/3, 2/5 <-> 3/5 and 3/7 <-» 4/7 
corresponding to flux 24, 41 and 58 0o respectively. Note that 
our new trial wavefunctions show an excellent reproduction of 
the increasingly nontrivial features as v — 1/2 is approached. 
In the top panel the data points have been shifted slightly 
horizontally from integer angular momentum values for bet- 
ter distinguishability. 



introduced here. 

The excellent agreement of our negative p wavefunc- 
tions with particle-hole conjugation of positve p wave- 
functions extends to the excited states, generated as out- 
lined in section |H] As examples, Figure ^ shows excita- 
tion spectra for v = 1/3, 2/5 and 3/7 and their respective 



particle-hole conjugate states. We also give, in Table |U 
values for magnetoroton gaps (which are the lowest en- 
ergy neutral excitations) as well as large k gaps (which 
are presumably the transport gap). These results, along 
with the excellent results for the ground-state energies 
confirm the validity of our approach to calculating CF 
wavefunctions at negative effective flux, which enables us 
to consider in the following section filling fractions above 
1/4 in the series v — p/{Ap +1), that were previously 
inaccessible. 



IV. RESULTS FOR 



: p/(4 P - 



The 4 CF series of composite fcrmion states, corre- 
sponding to filling fractions around v — 1/4, has been 
the subject of some recent experimental work|8j, yet the 
branch of filling fractions above 1/4 has been mostly in- 
accessible to numerical investigations. The numerical ap- 
proach we took for examining these states permits us to 
calculate the excitation spectra for systems with a mod- 
erate number of effective magnetic flux quanta N^ s = 2q 
for the CF-system. For negative p, the calculational com- 
plexity of the wavefunction increases with q, consequently 
the achievable system size is reduced compared to the sys- 
tems with positive p. For this very reason, extrapolation 
of the results to infinite system size in order to obtain the 
gap in the thermodynamic limit, and thus the mass of the 
CF, is difficult for these states. Nonetheless, the data ob- 
tained so far reveals that the spectra of states with p and 
— p at the same system size appear to be more and more 
similar as \p\ increases and v = 1/4 is approached. Table 

summarizes the numerical results for the groundstate 
energies and gaps of states around v = 1/4, and Figure 
121 shows several examples of excitation spectra revealing 
the above mentioned similarity. 

It is interesting to note that the p <-*■ — p similar- 
ity in spectra found here differs from the particle-hole 
symmetry correspondence of p with the — p — 1 states 
in the v — p/(2p + 1) series shown above in section 
IIIII One might naively expect that the particle-hole 
symmetry p <-» — p — 1 for states in the 2 CF series 
v = p/(2p + 1) above survives the composite fermion- 
ization attachment of two more Jastrow factors and re- 
sults in an approximate symmetry for v — p/{Ap + 1) 
states. For example, we could take the symmetry related 
m = 1 pair (p = -3,N = 12) and (p = 2,N = 10) 
at — 21 on Table |U and add two Jastrow factors 
so Ncf, —> Ncf, + 2(N — 1) resulting in an approximate 
symmetry relating (p = —3,N = 12, = 43) to 
(p = 2,N = 10, N<f, = 39) in Table HJ However, when 
we examine the overall shape of the resulting excitation 
spectra, we find very little relation between the two spec- 
tra related by p <-» — p — 1 for 4 CFs. In contrast we see in 
Figure [21 that the p <-> — p related states have quite sim- 
ilar dispersions, and the similarity appears to increase 
with increasing \p\. 

We comment that another interesting test of our ap- 
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-0.38266(3) 


0.0154(7) 


0.0169(7) 




2 


58 


-0.36765(3) 


0.0137(6) 


0.0137(6) 


18 




-3 


65 


-0.38126(3) 


0.0140(8) 


0.0190(8) 




3 


71 


-0.36427(3) 


0.0119(7) 


0.0171(6) 


21 




-4 


76 


-0.38017(3) 


0.0126(8) 


0.0192(10) 




4 


84 


-0.36182(2) 


0.0120(8) 


0.0175(8) 


24 




-5 


87 


-0.37940(3) 


0.0126(10) 


0.0196(9) 




5 


97 


-0.36006(4) 


0.012(1) 


0.017(1) 


27 




-6 


98 


-0.37883(3) 


0.0134(14) 


0.0195(14) 




6 


110 


-0.35869(2) 


0.010(1) 


0.016(1) 


oo 




-oo 


oo 


-0.37404(3) 


0.010(1) 


0.0238(8) 




oo 


oo 


—0.34842(4) 


0.0088(12) 


0.021(2) 


20 


-4 


-1 


75 


-0.37206(3) 


0.0123(6) 


0.0123(6) 


4 


1 


77 


-0.36788(3) 


0.0109(7) 


0.0109(7) 


24 




-2 


90 


-0.37172(4) 


0.010(1) 


0.0117(10) 




2 


94 


-0.36491(3) 


0.0118(9) 


0.013(1) 


28 




-3 


105 


-0.37165(4) 


0.010(1) 


0.0148(7) 




3 


111 


-0.36285(3) 


0.011(1) 


0.012(1) 


32 




-4 


120 


-0.37142(3) 


0.010(1) 


0.013(1) 




4 


128 


-0.36133(3) 


0.0093(9) 


0.014(1) 


36 




-5 


135 


-0.37138(2) 


0.012(1) 


0.016(1) 




5 


145 


-0.36021(3) 


0.011(1) 


0.014(2) 


oo 




-oo 


oo 


-0.3705(1) 


0.010(3) 


0.019(3) 




oo 


oo 


-0.35129(5) 


0.0087(19) 


0.017(2) 



TABLE II: Numerical results from Monte Carlo calculations of 4 CF wavefunctions for groundstate energies and gaps at filling 
fraction v = p/(4p +1), given in units of e 2 /e^o- Alongside the newly calculated negative p Quantum Hall states of this series 
above v — 1/4, we give the results for states with the respective positive p and equal particle number for comparison. We find 
that the similarity between the results for p and —p increases with increasing p. Two excitation energies are given from the 
spectra of neutral excitations: the magneto-roton gap as the lowest lying excitation and the large k gap measured at the highest 
possible angular momentum. Extrapolation to the thermodynamic limit has been performed using the above given data by a 
simple linear regression over the inverse particle number (see Figure [3J. In the case of groundstate energies, this extrapolation 
is based on density corrected values^J (Eg OTT = Eg^JuN^/N). 



proach can be obtained by constructing a trial wave- 
function for v = 1/3 as the p = — 1 member of the 
v = p/(4p + 1) series. We find that the ground state 
energy of this wavefunction is, within numerical preci- 
sion, precisely the same as that of the Laughlin v = 1/3 
trial wavefunction, which leads us to believe that we have 
precisely constructed that state. Similarly, we can exam- 
ine the excitation spectrum of v — 1/3 by using a 2 CF 
wavefunction of the series p/(2p + 1) with p = 1 or by 
using a 4 CF wavefunction of the series p/(ip + 1) with 
p = — 1. We find that the spectra obtained in these two 
approaches are quite similar to each other (albeit not 
quite identical), which gives us still further confidence in 
our approach. The values of the magneto-roton gap ap- 
pear to be almost exactly the same in both cases. Yet, 
examining the large k gap, the 4 CF trial wavefunction us- 
ing negative flux (p — —1) yields a slightly larger value, 
which decreases slightly more rapidly with the system 



size though, and seems to extrapolate to almost the same 
value at infinite N (see Figure [3}. 

Figure [3] shows the extrapolation to the thermody- 
namic limit of the data given in Table HT1 comparing gaps 
at p with those of states at ~p. A comparison of the ex- 
citations at v = 1/3 considered in the two different man- 
ners discussed above is also displayed there. Further, this 
figure gives us an idea of the quality of our extrapolation. 
Where the extrapolation is not particularly smooth, we 
cannot claim to deliver more than a rough result. The 
extrapolation of p = ±2 appears to be the most difficult, 
since it is not easy to distinguish a clear linear depen- 
dence of the large k gap as a function of N -1 for the 
initial data set. The reason for this problem appears to 
be that the magnetoroton gap is located at a large value 
of L, so that in the smaller systems, the magnetoroton 
gap is located very close to the largest k available, or even 
coincides with this point. For the p = 2 case, we can see 
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FIG. 3: This figure illustrates the extrapolation of the gaps 
to the thermodynamic limit for different filling fractions in 
the series p/(4p + 1) by means of a simple linear regression of 
the available data points over the reciprocal particle number. 
States with p and — p are displayed together. The state at 
v = 1/3 can be obtained in two different manners, i.e. as 4 CF 
in negative flux or as 2 CF in positive flux, which accounts for 
the two additional sets of data in the upper left. The mag- 
netoroton gaps of these different v — 1/3 states are close to 
indistinguishable, whereas the large k gap is slightly bigger for 
negative effective flux. As discussed in the text, extrapolation 
is least certain for p = ±2 due to finite size effects. 



FIG. 2: The method developed here permits us to calculate 
the dispersion curves for the low energy excitations at filling 
fractions p/(4p + 1) for negative p, i.e. above a quarter fill- 
ing, that were previously inaccessible. A comparison to the 
spectra of states below and above v = 1/4 shows, that the ex- 
citation spectra for p and — p show very similar features and 
seem to become more similar as p is increased. As above, the 
p > spectra are calculated using the method of Jain and 
Kamilla 4] whereas the p < spectra are calculated using the 
method discussed in the current paper. 



the error in this extrapolation clearly by comparing our 
extrapolated result to a result of a similar calculation us- 
ing larger system sizes from Ref. 0. This comparison is 
show in Figure 01 (see below). We might guess that the 
error in extrapolation for p — —2 is of similar magnitude. 

In order to obtain the composite fermion effective 
mass, we equate the activation gap A (determined from 
the excitation energy at the maximum angular momen- 
tum, i.e. biggest particle-hole separation) to the cy- 
clotron energy of CFs in their effective magnetic field: 



MB, 



ell' 



\2mp+ l\m* t\ 



(9) 



Since the gap is measured in units of the Coulomb in- 
teraction, we write A = (e 2 /e£ )S. Further taking 
into account e r — 12.8 for GaAs and the free electron 



mass as our point of reference, we find the dimensionlcss 
normalized 18] effective mass m* or = m* / (m e \fW v [T]) 
to be given by m* or = 0.0264/(|2mp + 1|<J). The 4 CF- 
masses we obtain are displayed in Figure 01 together with 
the large k gaps from which they are deduced. In addi- 
tion, for v < 1/4 (positive p) we have shown data from 
Ref. |a where larger systems were used than we have used 
here (When we also use larger system sizes, our results 
agree very will with those of Ref. |y). For v = 1/3, as dis- 
cussed above, we have shown results that treat this either 
as a p = — 1 point of the pj (4p+ 1) series or a p = 1 point 
of the p/(2p+ 1) series. For the latter case, we have like- 
wise included the result from the above mentioned larger 
system calculation. 

In Figure 01 we have intentionally displayed data ex- 
trapolated from a restricted set of small systems sizes 
with positive p so as to match the same set of system 
sizes that we study for negative p where we cannot go to 
very large systems. One can estimate the finite size error 
for our negative p calculations by examining the devia- 
tions between these restricted extrapolations at positive 
p compared to the larger system calculations of Ref. @ 
(also displayed in our figure). 

The behavior of the negative p fractions seems to 
roughly mirror the behavior of the positive p fractions. 
As was seen previously in Ref. for p > it is seen that 
the effective mass increases with \B e g\. As v — 1/4 is 
approached from either side, it is not clear if the effec- 
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// = 1/3. Below quarter filling, one observes a continuous 
increase of this mass. A similar trend appears at small 
values of the effective magnetic field above this point, 
but then the mass drops down again at v = 1/3, the 
final point of this series. 



Generally, our values for the effective mass seem to be 
lower than those measured in the experiments of Ref. 
by a factor of roughly 2.5. This error is rather expected, 
since similar discrepancies have been observed in previous 
studies based on the composite fermion picture 6] . It is 
known, however, that taking into account the finite width 
of the 2D electron gas changes the interaction so as to 
increase the effective mass [(J. 



FIG. 4: Masses of the composite fermions near v = 1/4 
and extrapolated large k gaps from which they were deduced 
at filling fractions of the series v = p/(4p + 1). Data from 
Ref. is from larger system size calculations and should 
be considered to be more accurate. Comparing these more 
accurate results to our smaller system calculations gives us 
an estimate of the finite size errors of our results which pre- 
sumably will hold even in the p < case. Except for the 
marked asymmetry between v = 1/3 and v = 1/5 we do 
not see much sign of an asymmetry around v = 1/4 which is 
observed experimentally 0. We also do not see signs of the 
experimentally observed diverging effective mass as v = 1/4 
is approached. Our values of the effective mass are somewhat 
less than that obseved experimentally in general. However, 
it is known that finite well width corrections tend to increase 
the effective massesp]. Two additional data points are shown 
in this figure (not mentioned in the legend) which give results 
obtained by treating v = 1/3 as a 2 CF state. 



tive mass will converge to a constant as would be pre- 
dicted by theory (up to logarithmic corrections pTjL In 
the experiments of Ref. Q, a striking asymmetry of the 
effective mass beween the high field and low field sides 
of v = 1/4 has been observed. While we cannot rule 
out some asymmetry from our data, we certainly cannot 
claim to see the extremely strong differences that are ob- 
served experimentally. This, however, is not surprising. 
Experimentally, the asymmetry is attributed to the prox- 
imity of a Wigner crystal state@- Since we are using a 
trial wavefunction approach, we should not see the effects 
of any imminent phase transition. 

Perhaps the most interesting data point in Figure 0] is 
the one at v = 1/3. Whether we treat this point as the 
p = — 1 member of the series p/(4p + 1) or the p = 1 
member of the series p/(2p+ 1), we find almost identical 
results of a very large gap, which establishes a continuity 
between the CF masses around v — 1/2 and v = 1/4. 
Furthermore, we note that this point is quite asymmetric 
with its reflection at v = 1/5. Certainly the hypothesis 
of constant effective composite fermion mass does not 
extend all the way from v — 1/4 out to both v = 1/5 and 



We note that we are not able to find any evidence of 
the divergence of m* as we approach v = 1/4 from ei- 
ther side, which is observed experimentally in Ref. 0. 
This is not surprising for several reasons. First of all, the 
experiment only sees strong divergences extremely close 
to v = 1/4 - which we cannot access numerically. Fur- 
thermore, one might suspect that the disorder might be 
the source of the divergences in the measurements used 
by Pan. More importantly, however, even if there were 
genuine infra-red divergences 01 of the effective mass as 
v = 1/4 is approached, one would not necessarily expect 
such divergences to be properly represented in a trial 
wavefunction approach. 



We emphasize that the most important achievement 
of this paper is not any particular numerical result. If 
genuine numbers were desired for comparison to exper- 
iment, we would want to use a more realistic interac- 
tion, accounting for finite well width[||, as well as per- 
haps Landau-level mixing, and we would want to use a 
more powerful computer to analyze ever larger systems. 
Instead we would like to emphasize in this paper that 
we have clearly demonstrated that we can extend the 
approach of Jain and KamillaQ to treat negative ef- 
fective magnetic field, and we can study these negative 
p composite fermion wavefunctions for reasonably large 
systems, which has not been done before. In this paper 
we have tested this method by using particle-hole symme- 
try for the 2 CF series and we have found our method to 
be quite accurate. We have then applied this method to 
4 CFs to study, for the first time, large systems for filling 
fractions p/(4p +1) with negative p. Our main physical 
result is that the effective mass appears to be roughly 
symmetric around and close to v = 1/4 although larger 
system calculations would be desirable. 

The authors acknowledge helpful conversations with E. 
H. Rezayi. G.M. acknowledges support from the French 
Ministry of Science, and thanks both the Ecole Doctorale 
de Paris and Lucent Technologies for their support that 
made participation in this research project possible. 
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APPENDIX A: CF WAVEFUNCTIONS WITH 
NEGATIVE EFFECTIVE FLUX 

The starting point for the composite fermion trial 
wavefunctions |3J are the single particle eigenfunctions of 
the quantum mechanical problem of a particle in a mag- 
netic monopole field on a sphere. (We use spherical co- 
ordinates with the azimuth 6 ranging from to 7r, and 
(j) the longitude ranging from to 2ir). The monopole 
harmonics are given bv|lOj 

K n %(fi) = 2 m M„(l - x) a ' 2 {l + xf/ 2 P^(x)e im * 
with a = —q — m, /? = q — m, g = \q\ + n + m, x = cos 



'2|g|+2n + l(|«| 



-m)l(\q\ 



m)\ 



4tt 



n\{2\q\ 



and Pg'P{x) are the Jacobi polynomials. This monopole 
harmonic represents an eigenstate of a particle on a 
sphere in a radial magnetic field with 2q flux quanta 
penetrating the sphere, where a positive sign refers to 
outwards pointing flux. Here, the angular momentum of 
the eigenstate is I = \q\ + n and the z component of the 
angular momentum m 6 {—I, . . .,1}. Further, n is the 
LL index n = 0,1,.... The above expression assumes 
the Haldane gauge [9j, where the singularities of the vec- 
tor potential are chosen to be located symmetrically on 
both north- and southpole of the sphere. Here, we fo- 
cus on the case where q < 0, since the a > case has 
already been discussed in detail in Ref. y. For the rest 
of this paragraph, we thus assume q < 0. Expanding the 
Jacobi polynomials in terms of the spinor coordinates 
u = cos(6'/2)e~ i ^ 2 and v = sin(6»/2)e^/ 2 , one obtains 

Y^(Sl) = (-l) n M q , n<m (u*)- q+m (v*T q - m x 

d- i )-(:)( i ^ i :;»)<»*»)-<«*»)"- 

This expression can equally well be obtained from 
the relation for complex conjugation of the monopole 
harmonics [TT| . if one corrects Kamilla's formula by re- 
placing q by \q\ in the appropriate places. We now use 
the Y q m as single particle wavefunctions (written as ipi 
in the main text) and composite fermionize by attach- 
ing Jastrow factors. As discussed in the main text, we 
can construct the many particle composite fermion trial 
wavefunctions by bringing the Jastrow factors inside the 
Slater determinant of single particle states, which may 
then be projected individually. Note, that in the spheri- 
cal geometry, the Jastrow factor becomes 



IjVk) 



(A2) 



The details of 
next. 



the required projection 
. . , un, vn)) are discussed 



First, we remark that the Jastrow factor^^ Jf is a 
LLL function, with q' = p(N — 1) zeros in Uj, i.e. it is a 
LLL wavefunction for flux q' > 0. Since iV is in general 
a big number, we have q' 3> \q\. The resulting wavefunc- 
tion, in turn, has to be a valid wavefunction for a total 
number of flux Q = q + q' > 0. Secondly, since projec- 
tion is a linear operation, we may consider the action of 
projection on each of the basis states Y^ m , separately, by 
expanding Jf in this basis. In general multiplication by a 
basis state Yfi followed by projection can be described 
as a linear operator called hereafter 2)^„. m . 



T>yi Y q = 91 9 Y q 

1 n,m 0,m' "r/q,n,m O.m' 



(A3) 



Since we know the entire basis of the subspace that we 
project upon, namely the LLL for flux Q with states \M), 
and \M\ < Q, the projection operator is J2 M \M}(M\. 
We now show, how this leads to an expression for 2) q 



q,n,m 



as a differential operator in the coordinate representa- 
tion, in which (IA3I) becomes 



£ Y %(n) dn>Y %*(n')Y* m (n>)Y q m ,(n>) 

M=-Q J 

= %i n , m Y( m ,{V) (A4) 

Integration over the longitudinal angle (f> singles out one 
nonzero scalar product for M = m + ml , and the one re- 
maining integral over the azimuthal angle 9 yields a well 
known binomial coefficient. Simplifying the normaliza- 
tion factors of Y^ m , on both sides, we have: 



n\ ( 2\q\ + n 



s J \\q\ + rn + s 



x(Nq, ,m) 4?r 



2 (q' - m' + s)\(q' + m' + n - s)\ 



(2q' + n + l)\ 

xuq -m+ q '-m' v q+m+ q '+m' = ^ q qnm u q '' m 'v q ' +m ' (A5) 



Using the explicit form of the normalization 
, 2 (2Q + 1)! 



M } 



Att(Q + M)\(Q - M)\ 



(A6) 



and remarking that the fractions of factorials that are left 
in this expression equal those that appear by the multiple 
derivation of a monomial u k 

Tu) ~ (Q + M)! U (A7) 

we may deduce 2)g n m by comparison of both sides: 
(2Q + 1)! 



21" 3 

*/ q ,n,m (2g/ + „ + !)] 



(-l)"M ? ,„, m X 

2\q\+n 



^ ^ ^3) [}q\ +m + sj * 



\q\+m+s 



\q\ — m+n — s 



v n - s (A8) 
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Let us remark that this result reproduces the known re- 
sult, that the projection on the LLL is achieved by per- 
forming the habitual procedure of moving all w*'s and 
v*'s to the far left, and replacing them with derivatives 
according to 



9 A 

— , and 

ou 



d_ 

dv 



(A9) 



Nevertheless, performing this explicit projection gives us 
a supplementary information in the form of a weight fac- 
for the different Landau levels before pro- 



tor 



(2Q+1) 



(2g'+n+l)! 

jection, which of course does not matter for the problems 
discussed here, but may play a role in other cases 15] . 

Practically, we would like to obtain a form of (|A8|) with 
the derivatives moved to the extreme right, which may 
be calculated using a straightforward application of the 
Leibniz rule for multiple derivatives in both u and v: 



E 



7 

(3 — a 



(A10) 



This yields a triple sum with the inner summation ranges 
being dependant on the outer summation index s. One 
finds that the summation ranges can be made indepen- 
dent of s since the summand becomes zero outside of the 
given intervals. As such the sum over s may be evaluated 
using 



E 

s=0 



— ( — l) a 8 n ,a+a' 



(All) 



Since the result yields a Kronecker delta, one of the re- 
maining sums becomes trivial, and after shifting the re- 
maining summation index, the final result is revealed to 
be exactly like (|A8|) . but with all derivatives placed at 
the very right. 

The projected composite fermion wavefunction is noth- 
ing but this operator applied to the single particle Jas- 
trow factor: 



Y q 

n.rn 



CF 



(A12) 



In order to perform numerical calculations with this 
wavefunction, we need to evaluate the derivatives explic- 
itly. One may use Jain and Kamilla's approach[3j, to 
commute the derivatives through the Jastrow factors as 



d_ 

dv. 



with 



&i = jr^jf, and Vi = jr^jf 

dm dvi 



(A13) 



(A14) 



Constructing a many-particle wavefunction out of these 
projected CF wavefunctions in the form of a Slater de- 
terminant, one may factor out the Jastrow factors again, 



and thus obtains a form which resembles single particle 
wavefunctions, on a basis of projected states Y£ with: 



(2Q + 1)! 
' (2q' + n 

E(-i) s f n V,f kl+n i 

^ \sj\\q\+m + s; 



(A15) 



Of course, this only appears to be a one particle wave- 
function, since there is an implicit dependence of the po- 
sitions of all other electrons in the system hidden in the 
operators Ui and Vj. The complexity of this expression 
increases with the total number of derivatives per term, 
given by Ng <0 — 2\q\ + n for negative q compared to 



q>0 _ 



n for positive effective flux. 



APPENDIX B: PARTICLE-HOLE 
CONJUGATION 

In the lowest Landau level, on the sphere, there are 
Nff, + 1 single particle eigenstates where N$ is the flux 
through the sphere. We label these eigenstates by the 
z-component of their angular momentum m. 

For a two-body interaction we can write the Hamilto- 
nian as 



H 



E 



Kti ,m 2 ,m 3 ,mi C mi c\ Jl2 C ma C m4 



(Bl) 



mi ,m2 ,m3 ,m4 



As usual, the normal ordering of the operators ac- 
counts for the uniform positive background. As usual, 
the fermion operators have anticommutation relations 
{c\,Cj} — Si t j. The interaction matrix V must have the 
following symmetries 



V, 



mi ,771-2 ,7713 ,7714 



= -V, 



ni2 ,777i ,7713,7714 



■V, 



= K 



7774,7773,7772,7711 



TTll ,7772 -7774 ,7773 

(B2) 



Furthermore, for any rotationally (translationally) invari- 
ant interaction, we must have angular momentum con- 
servation, which implies that the matrix element is zero 
unless 



mi + mi = nj3 + TO4 



(B3) 



We define the vacuum state |0 e ), to be the state which 
contains no electrons at all. The filled Landau level, 
which we can also think of as the vacuum for holes is 
written as \0h) = |Fillcd e ) = rim c m|0e)- It is convenient 
to introduce creation and annihilation operators d, for 
holes, given by 64 — cj , d\ = Ci which also obey the 
usual anticommutations {d],dj} = 5^. We now rewrite 
the Hamiltonian in terms of these hole operators. Us- 
ing the commutation relations as well as above described 
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symmetries of V we obtain 

H = 



^ . Kii : m2 ,m3 ,m4 ^mi ^L 3 ^m 4 



mi ,7712 i m 3 ? m 4 



2j2u m (l-2dld m ) + H d 



(B4) 



A general state containing N electrons in the LLL can 
be written as: 



■ |0 e ) (B9) 



where 



m,rri2 ,m2,m 



(B5) 



and 



/ / mi, 



m2 ,m3 ,m4 



mi ,m2,m3 ,7714 



We show below that for any rotationally invariant in- 
teraction, C/ m is actually independent of m. Furthermore, 
it is very easy to show that the energy of the entirely filled 
Landau level is given by 

^filled = 2 2^ frmi,m»,ma,TOi = (B7) 



Thus we have 



H = f 1 - HT^rM ^finod + H d (B8) 
\ iV0 + 1 J 

and JVft is the number of holes (i.e., the eigenvalue of 

Era dlndm). 



To particle-hole conjugate this state we construct 



I*) 



/ j "'mi . . . . ,771 iv 



4a 1 ---dl N \0 h ) 



(BIO) 



which is now a state containing ./V holes, or N$ + 1 — N 
electrons. Note that both of these states "live" in the 
same lowest Landau level which has + 1 single particle 
eigenstates (indexed by m). 

Now, if '5 is an eigenstate of H with eigenvalue 
(assumed to be real), then, since Hd has precisely the 
same structure as H we see that ^> is an eigenstate of Hd 
with the same eigenvalue. Thus, we obtain the particle 
hole conjugation relation 



E* = 1 - 



2N h 



Ued 



■Es 



(Bll) 



Lemma: U m is independent of m: 

If an interaction is rotationally invariant, we can choose 
any rotation R and write 



V mi ,m2,m 3 ,m4 — ^mi,mi (R) D m 2 ,m' 2 (R) Vm^m'^m'^m^ [D l m3 ^(R)}* [D l m4 m ^(R)}* (B12) 



where the -D's are rotation matrices as in Ref. Il2j and 
Z = 2^. Setting m2 = 777,3, and mi = 7714 = m and 
summing over 777.2 as prescribed in Eq. IB5l we obtain 



as well Eq. IB3I Since Eq. IB 131 must be true for any 
rotation, we can integrate over all rotations and use|12j 



U„ 



£ D km'S R ) Vm{,m- 2 ,m>s,m{ [D l m , m{ (R)}* f dR ^ ^ (R)[D l mm , (R)]* = constant (B15) 

n 1: m 2 J 1 1 



where we have used the orthogonality fl2^ 

^ Dm 2 ,m'y(R)[Dm 2 ,m' x (R)]* = ^ro!,,', 



(B13) 



(B14) 



independent of m, which shows that U m is independent 

of 777. 



[1] For a review of Composite Fermions see "Composite 
Fermions" , ed. O. Heinonen, World Scientific, 1998; and 
therein. 

[2] See J. K. Jain in "Perspectives in Quantum Hall Ef- 
fects" , ed. S. Das Sarma and A. Pinczuk, Wiley, 1997; 



and therein. 

[3] R. K. Kamilla and J. K. Jain, Phys. Rev. B 55, R4895- 
R4898 (1997). Int. J. Mod Phys. Bll, 2621 (1997). 

[4] R. K. Kamilla and J. K. Jain Phys. Rev. B 55, R13417- 
R13420 (1997). 



11 



[5] R. K. Kamilla, J. K. Jain, and S. M. Girvin Phys. Rev. 
B 56, 12411-12416 (1997); K. Park and J. K. Jain, Phys. 
Rev. Lett. 80, 4237-4240 (1998); K. Park, V. Melik- 
Alaverdian, N. E. Bonesteel, and J. K. Jain Phys. Rev. 
B 58, R10167-R10170 (1998); K. Park and J. K. Jain 
Phys. Rev. Lett. 83, 5543-5546 (1999); V. W. Scarola, 
K. Park, and J. K. Jain Phys. Rev. B 61, 13064-13072 

(2000) ; K. Park and J. K. Jain Phys. Rev. Lett. 84, 5576- 
5579 (2000); T. Sbeouelji, K. Park, J. K. Jain, and N. 
Meskini Phys. Rev. B 62, R4802-R4804 (2000); K. Park 
and J. K. Jain Phys. Rev. B62, R13274-R13277 (2000); 
V. W. Scarola, K. Park, and J. K. Jain Phys. Rev. B 
62, R16259-R16262 (2000); S. S. Mandel and J. K. Jain, 
Phys. Rev. B 63, 201310 (2001); Phys. Rev. B 64, 081302 

(2001) ; Phys. Rev. B 64, 125310 (2001); Phys. Rev. Lett. 
89, 096801 (2002); T. Sbeouelji and N. Meskini, Phys. 
Rev. B 64 193305 (2001). 

[6] X. Zu, K. Park, and J. K. Jain Phys. Rev. B 61, R7850- 
R7853 (2000). 

[7] K. Park, N. Meskini, and J. K. Jain, J. Phys. Cond Matt 
11, 7283 (1999); V. W. Scarola, S. Y. Lee, and J. K. Jain, 
Phys. Rev. B66, 155320 (2002); K. Park and J. K. Jain 
Phys. Rev. Lett. 81, 4200-4203 (1998). 

[8] W. Pan, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, K. W. 
Baldwin, and K. W. West Phys. Rev. B 61, R5101-R5104 
(2000). 

[9] F. D. M. Haldane, Phys. Rev. Lett. 51, 1583 (1983). 

[10] T. T. Wu and C. N. Yang, Nucl. Phys. B107, 365 (1976). 

[11] T. T. Wu and C. N. Yang, Phys. Rev. D 16, 1018 (1977). 

[12] A. R. Edmonds, Angular Momentum in Quantum Me- 
chanics, 2ed, (Princeton Press, Princeton, N.J., 1974). 

[13] R. Morf, N. d'Ambrumenil, and B. I. Halperin, Phys. 
Rev. B 34, 3037 (1986). 

[14] B. I. Halperin, P. A. Lee, and N. Read, Phys. Rev. B 47, 
7312 (1993). 



[15] G. Moller, S. H. Simon, to be published 

[16] We note in passing, that one could certainly construct 
several other trial wavefunctions which project in slightly 
different ways. For example, for m > 1 we might consider 

^(fj) = ? {i'i('> f j)Jj 11 } witl1 mi + m 2 = m (and 
mi > 1). In the cases we have checked we have found 
that wavefunctions built with this version of ip are also 
extremely similar to those built from Eq.HJ 

[17] A subtlety in this process results from the observation 
that the wavefunction, calculated as described in the ap- 
pendix [X] by separating it into a Slater-determinant of 
the projected pseudo single-particle-wavefunctions and a 
Jastrow-factor, tends to produce numerical instabilities 
in the determinant algorithm. This is easily understood, 
since two particles approaching each other at a distance 
d will cause the respective to grow as oc op(l' ? l+ n > 
in the same manner, producing two linearly dependant 
columns in the matrix. Now, due to numerical errors, 
instead of obtaining zero when two columns become the 
same, the evaluation yields numerical errors following ap- 
proximately the same power law behavior, that will dom- 
inate the result. Since the short distance behavior of the 
wavefunction is known, extrapolation allows us to correct 
for such events that occur with a probability of less than 
10~ 4 in all of the presented calculations. 

[18] Note that Pan et. al. 8] use a a different normalization 
for m* , which amounts to a slight rescaling as follows: 

m Pan = ^jBv/Bx/imlar. 

[19] in this appendix, we note the number of flux pairs at- 
tached to each CF as p, according to Kamilla's notation, 
whereas m, that had been used for this purpose in sec- 
tion[H]to comply with the common notation for the filling 
fractions of Jain's series, is already used for the eigenvalue 
of L z here. 



